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Abstract 

We calculate renormalization group flow equations for the linear a- 
model in large Nc approximation. The flow equations decouple and can 
be solved analytically. The solution is equal to a self consistent solution 
of the NJL model in the same approximation, which shows that flow 
equations are a promising method to extend the calculation to higher 
order in l/Nc- Including explicit chiral symmetry breaking, the large 
Nc approximation describes physics reasonably well. We further compare 
the analytic solution to the usually used polynomial truncation and flnd 
consistency. 
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1 RG flow equations of the linear cr-model 



The physics of strong interactions in the nonperturbative strong couphng regime 
cannot be obtained directly from QCD. Especially, it is not yet possible to 
integrate out the gluonic degrees of freedom in a rigorous way Q. From 1- 
gluon exchange or instanton models one is lead to a local 4-fermion interaction 
described by the Nambu— Jona-Lasinio (NJL) model. The NJL model is non- 
renormalizable and therefore includes an explicit cutoff scale. It is commonly 
considered in large Nc approximation . If fluctuations in the quark condensate 
are taken into account one needs an additional cutoff which reduces the 
predictive power of the model. Such an additional cutoff is not needed in the 
renormalization group (RG) framework, when the NJL-model is embedded into 
the linear cr-model. The main aim of this work is to clarify the connection 
between the NJL model and the RG flow equations for the associated linear 
cr-model in the large N^. approximation. 

RG flow equations describe the average of an effective action and represent the 
continuum analogue of a block spin transformation. Various kinds of exact flow 
equations have proved successful in describing a huge set of different systems 
[|j . For an introduction to the RG method we refer to the review [|j . The flow 
equations are ultraviolet (UV) and infrared (IR) finite through the introduction 
of a scale (fc) - dependent cutoff function. Degrees of freedom with momenta 
larger than the coarse-graining scale k are integrated out. The solution of 
the flow equations includes summations of different diagrams and is similar 
to the solution of a coupled set of Schwinger Dyson equations. We use the 
Schwinger proper time regularization method Q , which preserves all symmetries 
of the theory and keeps the physical interpretation of the evolution equations 
particularly simple for phenomenological applications. This simplicity makes 
these RG improved 1-loop equations interesting, although they are not exact 
RG flow equations 0. 

The lagrangian of the NJL model with an explicit current quark mass rric reads 
in Euclidian spacetimc 

C^,^ = + mc)q - G {{qqf + {qif-fsqf) . (1) 

A partial bosonization by the introduction of auxiliary scalar fields cr and tt 
gives the associated linear a-model 

= 0q + gq{a + itttj^) q + — {a^ + Tr"^) - 6 a . (2) 
The respective couplings are connected by 

2G = ^ , TOc = ^ . (3) 

In order to solve the theory it is sufficient to compute the effective action, which 
is the generating functional of one particle irreducible Feynman diagrams. It 
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involves in principle arbitrary higher order terms consistent with the symmetry 
of the initial action. Using RG methods we are going to calculate the effective 
action in a suitable truncation. Surely we cannot account for all possible terms, 
so our approximation is to assume that the linear cr-model with quarks is a valid 
description of nature at a scale of 0.2 GeV<A: < 1 GeV. In this region the gluon 
degrees of freedom are supposed to be already frozen out, but confinement is 
not yet relevant. It remains to be investigated whether the effects of other qq- 
resonances are sufficiently taken into account by quarks coupling to (a, 7?) fields. 
The physics of the model is governed by chiral symmetry, which is spontaneously 
broken in the IR. The Euclidian UV action of the general linear cr-model is given 

by 



^[$, g, = j d^x (^Zg 0q + +9qMq + U 



(4) 



The quark- and meson fields g, M are bare fields, where $ — {a, vr) is the 
0(4)- and M = (j + ifn'y^ the chiral SU{2)l 5'[/(2)/j-representation of the 
meson fields. We consider this action to be preserved during the evolution and 
choose the starting values 

Z^ = 0, (5) 
Zg = l, (6) 



and the UV potential 



™2 

t/uv-^1>'-5a, (7) 



in order to obtain Eq. (^) as UV lagrangian. For comparison with the standard 
NJL approximation we will work in the following in the chiral limit riic ^ S = 
0, but since the term linear in a does not get renormalized Q, the following 
computation is also valid for explicit symmetry breaking. This case will be 
discussed in Sec. 3. The effective action F, which depends only on averaged 
fields, can be computed by doing a saddle point expansion corresponding to a 
one loop approximation (cf. [^): 

1^ . /(5 ^S[<i>,q,q] \ 



F[<i>,,,,].5[<I>,,,,]--Trlog^^^^j (8) 

^ ^^S[<i>,q,q] / ^^S[<l>,q,q] y^ J^S[<i>, q, q] \ 
(5$i(x)<5$j(y) 5$'(x)<5q(y) UqW<5q(y)y <5q(x)5$j(y) j ' 

The first logarithm results from fermion loop fluctuations, whereas the two terms 
in the second logarithm are the contributions of the bosonic and the mixed loop 
respectively. 
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The logarithms are transformed into Schwinger proper time integrals regularized 
by a heat kernel cutoflF /(fc^r) (cf. |^) 

log(A) = -y^ —e-^H{]<?T). (9) 
As cutoff the following set of functions is used 

/(")(fcV)^X^&l^exp(-fcV). (10) 

The heat kernel cutoff functions Z^"-* suppress fluctuations with momenta below 
the cutoff scale k. Going to fc ^ means to include more and more IR modes. 
Finally all modes are included, because of the limiting behaviour /'"^ — + 1. Due 
to the extra r terms in front of the exponential the integral diverges in the UV 
region only through an additional term ~ log A^ . To keep the following compu- 
tation simple, we will consider the cutoff /'^^•' used in It has been shown for 
the case of a 0(A^)-model, that the RG equation derived with the cutoff /^^^ is 
equivalent to an exact RG equation ]lO| . 

The detailed computation of the effective action is based on a derivative expan- 
sion and will be presented in a separate paper The purely bosonic part 



was investigated in |12 



The derivative of the regularized effective action with respect to the cutoff scale 
k yields a flow equation for the effective action. Identifying terms in ^ with 
similar terms in S'[$, g, q] we obtain evolution equations for the effective poten- 
tial [/, coupling constant g and the wave function renormalization factors Zq 
and 

To improve the one loop equation we substitute the masses and couplings of 
the classical action by the running masses and running couplings of the effective 
action T. This replacement turns the one loop equations into RG improved flow 
equations, which include higher loop terms successively into the proper time 
integral. Therefore the flow equations can treat nonperturbative physics in the 
strong coupling region [B . 
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We find the following coupled set of equations for the parameters of the model|^ 

ac/ ^ _ iv^ j_ ( k^ k^ \ 

dk 87r2 k-^ + M^ 32Tr^ \k^ + My k^+M^J' ^ ' 

Bg 1 .g3 / fc6(2fc2 + M2 + M2) 

k-r- — ■ ' 



dk 167r2 Z^Z^ Uk'^ + M^)^{k^ + M^ 



2\2 



{k^ + M^)^{k^ + M^y 



(12) 



k^ 

" dk 




fc6 (2fc2 + M2 + M2) 



(fc2+Af2)2(fc2 + M2)2 



(13) 



dk 47r2 Z2 (fc2 + M2)3 87r2 * (k^ + M^ik^ + M^)^' ^ ' 

where the effective masses are given by 

^,2(^2. A ^ 

' Z$ ,5$2 ' 
and the effective 4-meson coupling reads 



2 ,5^[/(^ 2) 



A('J'')-^^77^- (16) 



In this set of equations all the masses and couplings are functions of the scale k. 
The heat kernel method gives these equations in a very simple form. Since the 
IR parameter k acts like a mass cutoff, the propagator structure of the different 
loop contributions is still recognizable. 

Although the effective masses and the effective 4-meson coupling in Eqs. ( plf - 
are fully renormalized quantities, the potential U and the Yukawa coupling 

^ Contrary to the equations for C/(<I>^ ) and Z$ ) , which are valid for general <l?^-dependent 
functions, the given equations for the Yukawa coupling and the quark wave function renor- 
malization factor are only valid for ^-independent g and Zq. The general expressions get 
correction terms which do not contribute in large Nc- 
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g are still bare ones. In order to obtain physical couplings one has to express 
the theory by renormalized fields qr, qr and ^r. This is done by rescaling the 
fields to absorb the wave function renormalization factors in front of the kinetic 
terms: 

= $ , qj, = zjq. (17) 

where $0 denotes the vaciuim expectation value of the field $. The renor- 
malized potential is a function of the renormalized meson field Ur{^^) = 
U {{Z^{^q))~^ and the renormalized Yukawa coupling is related to the bare 
coupling as gR = (Z$($o))~2 g. Finally, the physical masses and the phys- 
ical 4-meson coupling are given at the vacuum expectation value of the meson 
field 

mi = M,($2) , ie{a,Tr,q} ; A = A($2) . (18) 
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2 RG in large Nc and NJL gap-equation 



In the large N^. limit only the fcrmion loop is relevant. Since the NJL effective 
action involves a trace over quark fields it is proportional to N^. and therefore the 
NJL coupling G is proportional to l/Nc- This also has to hold in the bosonized 
case Eq. (^). Assuming that the meson masses are independent of one 
obtains from Eq. (^) that g behaves oc 1/y/Nc, in the same way as the QCD 
coupling constant. This implies that the meson fields <& scale as 0{^/W^), which 
in turn means that the effective 4-meson coupling A should be proportional to 

Due to the large Nc behaviour, the quark contribution in the flow equation of 
the potential U Eq. ( pT|) is of order Nc, whereas the meson loop term is of 
order one and can be neglected. Also, for Z$ Eq. ( p^ ) the second term can be 
dropped since it is of order compared to the first one which is of order 

one. The wave function renormalization parameter Zg Eq. ( [l^ ) and the bare 
Yukawa-coupling g Eq. (|l2|) do not evolve, since the right hand sides of the 
evolution equations are of higher order in 1/Nc than the left hand sides, which 
reflects that the corresponding fluctuations are given by mixed loops. 
The large Nc flow equations read {Zq=g = l) 

dU ^ ^N^ 



dk 87r2 k'^ + M, 



dZ^ NfNc k^ 

dk 47r2 (fc2 + Af2)3 ■ ^ ^) 

The right hand sides of these equations only depend on the quark mass Eq. 
(^5|), but not on the potential U or the wave function renormalization parameter 
therefore they completely decouple and can be integrated analytically. In 
particular, these equations do not contain a RG improvement. Note that in the 
large Nc limit we have = <I>2 and the integrals read 

U{k) k 

dU = -^ N^d-^ (21) 

Z^ik) k ^ 

dZ^ = -^^^ / TT^-r^dn . (22) 
47r2 J (k2 + $2)3 ^ ' 

Neglecting ^-independent terms, the integral for the potential U yields 

L.(<i>,fc) = [/..($) + ^ [^'lo^{^^^)~^\kt-k^)) , (23) 

where the UV potential is given by Eq. (|^). Using = as starting value, 
the equation for Z$ can be integrated in the same way giving a slightly longer 
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expression 



87r2 V ^\k^ + $2 

(fcg, ~ k^) (4 e fcg, + 3 (fc^ + fcg,) $2 + 2 j>-^ 

(fc2 + $2)2 (^2^ ^ ^2)2 



(24) 



Eqs. ( |23| , 24) are the exact analytic solution of the large Nc flow equations for 
a given scale k. To obtain the physical vacuum expectation value <&o(fc) of the 
field, the potential has to be minimized with respect to <&: 



JC/($,A:) 



(5$2 



(25) 



This equation cannot be solved analytically due to the non algebraic logarithmic 
term. As we will show, it corresponds to the self consistent NJL gap equation 
at the IR scale k = 0, determining the physical quark mass niq = $o(A; = 0) 
as the minimum of the potential Eq. (|2^). The full potential provides more 
information than the NJL gap equation, since it gives arbitrary higher order 
moments of $ besides its vacuum expectation value. This allows to compute 
for instance the spectrum of the Dirac operator which will be presented in a 
forthcoming paper Q. 

The NJL model in large Nc approximation is standardly considered in a self 
consistent approach which leads to the NJL gap equation [||. This equation 
determines the constituent quark mass and is usually regularized by a hard 
cutoff. Furthermore, to make contact with chiral degrees of freedom there is an 
additional equation that determines the mesonic wave function renormalization 
factor once the constituent quark mass is known 

To compare with our flow equations we regularize these equations in a Schwinger 
proper time scheme using the heat kernel cutoff. Since the heat kernel cutoff 
/^"^(fc^r) acts like a IR cutoff function, the desired UV regulation is given by 
the function (1 — f'^^\k^^T)), which cuts off momenta larger than fcuv- The gap 
equation and the determining equation for Z$ read 

poo f 74 

1 ^ ^GNcNf dr e-(^'^+"?)(l - &Hklr)) , (26) 

^ 2NcNf rdr J ^ e'^(p''+<\l - &\kl^T)) . (27) 

The integrals appearing in these equations can be integrated analytically. With 
Eq. (|^) the NJL coupling G is connected to the UV meson mass m^^. For the 
cutoff /(^^ the integrated form of the gap equation Eq. (|2^) is analytically equal 
to Eq. (pS) . Similarly, the determining equation for Z$ Eq. (p7|) yields Eq. (e3) 



8 



evaluated at $^ = after integration. This connection between the equations 
in the RG scheme and the equations in the self consistent NJL approach holds 
for all cutoff functions Eq. (^0|), provided both schemes are regularized 
with the same cutoff. This shows that the RG flow result is exactly identical to 
the standard NJL approach in large Nc approximation. 

Although the vacuum has to be computed numerically, there are two quantities, 
that allow an analytic solution. The first one is the scale fc^ss, where the 
system evolving from the UV enters the chiral broken regime, which is given as 
the solution of the equation 



^$2 



= , (28) 



and takes the simple form 



kxSB = W - ]^7^"^uv • (29) 

The other one is the ratio ma-/mq, which results from Eqs. (|l^, |l8|) using the 
analytic expressions for U and Z$ at fc = and $ = $o- It is exactly 2 in large 
Nc, as known from the standard NJL approach 

For a numerical solution, we adjust the UV parameters to reproduce physical 
values of niq = 320 MeV and = 93 MeV in the case of an explicit chiral 
symmetry breaking, discussed in the next chapter. The initial parameters for 
the RG evolution read 

k^^ = 1037 MeV , 

= 228.0 MeV , (30) 

and Zq = g = 1 = const. This corresponds to a four fermion coupling of 
G = 9.62 GeV^^. With these initial values we calculate the effective potential 
from Eq. ( p3|) and the wave function renormalization from the corresponding 
equation for The result is the renormalized effective potential Ufj, which 
is shown in Fig. |l| as a function of the field <i>fl and the scale k. From this 
potential we compute the scale dependent pion decay constant f-r^ik) shown as 
white curve in Fig. ^ 

The quark mass and meson masses, which are given from Eq. ([T^ ) are plotted 
in Fig. H (solid lines). In the chiral symmetric regime for scales fc>fc^5B = 626 
MeV (cf. Eq. (|2g|)) both meson masses are identical, while the quark mass 
is zero. Below the symmetry breaking scale the pions become massless while 
the CT-meson mass rises again. Due to the nontrivial vacuum of the effective 
potential, the quarks acquire a finite mass. In the IR limit fc = we find the 
physical values for the the quark mass = 'I>o(fc = 0) and the u- meson mass 
Eqs. (p|,r 
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0.05 0.1 0.15 

[GeV] 

Figure 1: The evolution for the renormalized meson potential in large Nc 
approximation. In the UV the renormalized potential is infinite for all ^r due to 
the vanishing wave function renormalization constant. The quadratic potential, 
which arises when gets finite, evolves to a potential with broken symmetry 
in the IR. The white curve shows the evolution of the vacuum expectation value, 
which corresponds to — MeV at fc = 0, denoted by the endpoint. 

— 1/2 

The Yukawa coupling is given hy gR = Z^ (fc = 0) and the pion decay constant 

is /7r = ^y^(fc = 0) $o(fc = 0). The IR values are summarized in Table ^ (second 
column) . 

3 Explicit symmetry breaking 

Although the physical number of three colors is not really large, we know from 
standard NJL computations, that the large Nc limit gives a reasonable approxi- 
mation. We next compare our model with the physical case of an explicit chiral 
symmetry breaking due to a finite current quark mass rric- 
The value of the mean current quark mass rric — — 7 MeV is taken from 

chiral perturbation theory which has been performed at our UV scale of 1 
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Figure 2: The renormalized masses of the exact solution as functions of the 
scale k in the chiral hmit (sohd) and with an exphcit symmetry breaking due 
to a current quark mass of 7 MeV (dashed). The quark mass in the symmetric 
regime and the pion mass in the broken regime vanish in the chiral limit, whereas 
the pion acquires a mass of m-^ = llQ MeV in the case of a finite current quark 
mass. 



GeV, ahhough there is some uncertainty because of the different schemes. From 
Eq. (H) the parameter 5 can be computed resulting in 5 = 3.64- 10~'*(GeV)'^. As 
already discussed, the solution Eq. (|2^) obtains in the case of explicit symmetry 
breaking, with a term linear in a in the UV potential Eq. @). The vacuum ex- 
pectation value is now given by minimization with respect to a. The evolution 
of the renormalized physical masses is shown as dashed curves in Fig. |^, which 
exhibits a very reasonable splitting of the mass scales. There is a continuous 
transition in the k evolution from the current quark to the constituent quark 
regime, since the system is in a nontrivial vacuum right from the beginning of 
the evolution. The quark mass assumes its current quark mass in the UV and 
the pion looses its role as an exact Goldstone boson and acquires a small finite 
mass in the IR. 

The IR values of the physical quantities are shown in the first column of Table 
, The effect of the explicit chiral symmetry breaking is a minor change in the 
other parameters and a finite pion mass of 116 MeV which is reasonable close 
to the physical value. The change of the pion decay constant compared to the 
chiral limit amounts only to 1 MeV and is rather small in large Nc approxima- 
tion compared to chiral perturbation theory computations |l4[. However, this 
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mc = 7MeV 


X-limit 


$'*-truncation 




[MeV] 


93.0 


92.0 


105 


A 




23.6 


22.6 


14.9 






3.44 


3.36 


2.73 




[MeV] 


320 


309 


285 




[MeV] 


650 


618 


570 




[MeV] 


116 









Table 1: The IR values for different approximations are shown. All computa- 
tions are done with the UV initial conditions kuv = 1037 MeV, ruuv = 228.0 MeV 
and Z^^^Q. The first column shows the analytic RG result with explicit chiral 
symmetry breaking due to a finite current quark mass of 7 MeV. In the second 
column the analytic RG result in the chiral limit is given, which is equal to the 
result of the NJL gap-equation regularized with the proper time regulator. In 
the last column the results of a simple quartic truncation of the potential are 
shown. 



seems to change at order 1/Nc pT[ ]. 

4 Truncated potential 

The analytic solution of Eqs. (|l^,^^ can test the usual parametrizations of the 
model. One relies on such parametrizations, if higher corrections are regarded, 
e.g. mesons are taken into account |^. A standard approach is to expand the 
model in powers of <I>^ up to a given order Pmax- 

The potential U evolves under the renormalization flow from a trivial to a 
nontrivial form. Therefore we consider two parametrizations. In the simplest 
quartic truncation with Pmax = 2, they read in the chiral symmetric regime 
(7($2^fc) = Ililhlq>2+Mk) ^^2^2 ti^g b^Qj^gj^ regime C/($2,/c) = ^($2_ 

$Q(fc))^. Since the wave function renormalization factor Z$ arises together with 
the dynamical term of mass dimension four, it is correspondingly expanded up 
to order Pmax— 2 and is therefore ^-independent in the simplest approximation. 
The evolution equations for the couplings and A can be derived from Eq. 
(^9|). We find for large Nc in the symmetric case {g = Zq = l): 

dm,^ _ NcNf 

'~dk' ~ 

d\ _ NcNf 

'dk ~ 2^ 



fc-T7^ = +^fc^ (31) 
= > (32) 
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and for the broken regime: 



d^l N^Nf 1 P ^^^^ 



(35) 
(36) 



dk 


47r2 


A(l 






NrNf 




1 


"dk " 


27r2 


(1 + 


^l/k^f ' 




NcNf 




1 


" dk " 


47r2 


(1 + 





To solve these equations numerically, we use again the initial values Eq. (po|). 
In addition we choose Auv = to have the same UV potential as before (cf . Eq. 

In the symmetric regime the exact and approximate solutions are identical. 
In the broken regime we switch to the second parametrization and follow the 
flow given by Eqs. ( ^4[]3^ ) to fc = 0. Now we have a set of coupled equations. 
This comes from the fact that we are tracking the evolution of the minimum $0 
of the potential U, which has no analytical solution from Eq. (p3|). In the IR 
we find a pion decay constant = 105 MeV compared to a value of 92 MeV 
for the exact solution (cf. Fig. ^). The IR masses and couplings from these 
truncated flow equations are given in the third column of Table ^ 



0.2 
. 175 
. 15 
. 125 
[GeV] 0.1 
. 075 
0.05 
0.025 



0.2 0.4 0.6 0.8 1 



k [GeV] 

Figure 3: The renormalized meson mass rriR = rria = in the symmetric regime 
and the pion decay constant in the broken regime as functions of k for the exact 
solution (solid) and the quartic truncation (dotted). In the symmetric regime 
the two solutions are equal, but the vacuum results for /j^ differ by 14%. 
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When the truncation is extended to higher orders Pmax, the IR value of /,r 
converges fast to the exact solution with a pion decay constant of 92 MeV, as 
indicated by the thick points in Fig. |^. On the other hand, if only the effective 
potential is expanded, the flow still converges, however to a strongly deviating 
result of 111 MeV, shown by the thin points in Fig. ^ This is in contrast to a 
naive dimensional analysis, which would suggest that the higher order momenta 
of Z$ should be suppressed by their high mass dimension in the vicinity of an IR 
fixed point. Although this does not seem to be the case in large N,, it may be the 
case when meson corrections are included, where a partial fixed point behaviour 
occurs. Nevertheless, in order to obtain a credible solution it is necessary to 
choose the truncation order in such away, that convergence with respect to all 
model parameters is achieved. Indeed, at order pmax = 7 the results are to the 
given accuracy exactly the ones of the exact solution in the second column of 
Table |l|. 
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Figure 4: The dependence of the physical pion decay constant /^r and the 4- 
meson coupling A on the truncation order Pmax of a polynomial parametrization. 
Thin points give the behaviour when only the potential U is expanded and the 
wave function renormalization factor Z$ is considered ^-independent. They do 
not converge to the exact results given by the dashed lines. Thick points show 
the dependence when U and Z$ are expanded and exhibit a fast convergence 
to the exact results. 
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5 Conclusion 



We have derived RG flow equations for the hnear cr-model in large approx- 
imation as a different approach to solve the NJL model. It was shown that in 
large Nc the RG approach reproduces the NJL gap equation and therefore in 
this limit the solutions of the two schemes are equal. In addition one finds from 
the flow equations the whole effective potential, instead of just its minimum. 
In the case of explicit symmetry breaking due to a finite current quark mass 
we find very reasonable meson- and constituent quark masses. A comparison of 
the full flow equations with a polynomial truncation has shown a deviation of 
« 15% in lowest order, but a fast convergence to the full result. 
Both, the exact and truncated calculation use two parameters fcuv and TOuv to 
fix the quark mass and the pion decay constant This is the conven- 
tional method in NJL theory. An advantage of RG equations is the reduction 
of the relevant UV parameters, if the theory exhibits a partial IR fixed point 
behaviour. Such a fixed point behaviour indeed occurs when meson loops are 
considered. 

In addition, meson loop corrections can be included systematically without the 
problem of introducing additional cutoffs, since the effective action is renormal- 
izable. The analytic large Nc solution presented here serves as a reference to 
the analysis of the full flow equations with meson fluctuations, which will be 
published in a forthcoming paper pT[ |. 
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